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I. INTRODUCTION 



Quantum fluctuations of surfaces are important in many physical phenomena. It is less 
known these fluctuations can sometimes be expressed in term of a functional integral over 
reparametrizations of variables relevant in the Feynman path integral. In this paper we shall 
consider the expression 

w(c) = Jve e ~ KA ^\ (1) 

where x(9) is the boundary curve C and where 

A[xm = jl r ia r d y m -m?, (2) 

1 V n 8tt J J 1 - cos(<7 - a') V ; 

The functional (2) is known in mathematics as the Douglas integral [1], whose minimum 
with respect to variations of the reparametrizations 9 gives the minimal area 

A[x(6.)] = S nin (C). (3) 

Here 9*(a) is the saddle point of the integral (1). The functional integration (1) thus gives 
the area law to leading order. 

The path integral over reparametrizations (1) was introduced in this context in Ref. [2] in 
connection with an off-shell string propagator. More recently it was proposed by Polyakov 
[3] as an ansatz for the Wilson loop in large N QCD. The leading behavior obviously gives 
the leading behavior of the Wilson loop, found in most string models, where the bulk field 
X m (t, a) satisfies the Dirichlet boundary condition 

^(r^)|boundary = ^((7). (4) 

To derive Eq. (3) from the Nambu-Goto action one can follow Douglas [1], or the more 
recent elegant paper by Migdal [4]. 

The functional integral (1) can be expanded around the saddle point 6*. With 9(a) = 
0*(cr) + (5(a) we obtain the first non- vanishing contribution 

A 2 = £ f\9 f\9> X (9). X (9>) My-M? (5) 
8tt J J yj 1 > i - cos (a*(9) - a*(9>)) yj 

Here a*(9) is the inverse function of 9*(a). 

Some dynamical consequences of Eq. (5) have been discussed by Rychkov [5] and by the 

authors [6]. However, it is perhaps fair to say that the physical meaning of the fluctuation 

integral (5) is not so clear. 



In this paper we shall show that the leading part of the fluctuations of large loops from 
(5) are transverse fluctuations of the minimal surface embedded by the curve x{6*). This 
we have shown for a rectangle and an ellipse, but we suspect the result to be more general. 

To explain our result we mention that in the Nambu-Goto action transverse fluctuations 
add a contribution 



(6) are called the Liischer term [7]. By lattice Monte Carlo calculations in three and four 
dimensions, this term has been found to occur [8] in quenched SU (N) for various N's. Thus, 
for large distances the two leading terms from the Nambu-Goto action describe QCD quite 
well. Therefore, for a rectangular boundary curve the T/R term can be identified with 
transverse fluctuations of the minimal surface. 

With these remarks in mind we now give our main result: after integration over the 
fluctuations (3 in the functional integral 



we obtain as the leading contribution the Liischer term corresponding to d=2Q. In general, 
we obtain the rf-dimensional contribution (plus the area term) from 



Here O is the operator which emerges in the semiclassical expansion of A [x(0)\ to quadratic 
order as exhibited in Eq. (5). This modification of Eq. (1) does not effect the classical 
limit leading to the area law, and has the meaning of a pre-exponential in the semiclassical 
approximation. Thus, this equation gives the leading effective QCD string behavior in terms 
of functional integration over reparametrizations. 

We mention that our considerations may have potential applications in condensed matter. 
For example, in the inverse square XY model one encounters expressions somewhat similar 
to Eq. (5), see for example Ref. [9]. We shall, however, not pursue this track in the present 
paper. 

The plan of this paper is the following: in Sect. II we discuss the general framework for 
the semiclassical approximation, and in Sect. Ill we carry out the functional integral over 
reparametrizations for a rectangle. A similar calculation for an ellipse is done in Sect. IV. 





(7) 




(8) 
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In Sect. V we derive the generalization given by Eq. (8) and in Sect. VI we make some 
conclusions. A number of more technical points have been discussed in some Appendices. 



II. A SEMICLASSICAL CORRECTION 

The path integral over reparametrizations (1) reproduces the exponential of the minimal 
area in the classical limit KS m i n — > oo. To calculate the semiclassical correction, we expand 

6(a) = 6.(a) + P(a) (9) 

or 

a(9) = a.(0) + 0(9) , (10) 

where 

0(0) = (3(2n) = (11) 

and expand the Douglas integral to quadratic order in (3 around the classical trajectory 
6* (a). This expansion makes sense because typical trajectories in the path integral over 
reparametrizations (1) are smooth as KS min — > oo and have Hausdorff dimension one [10]. 

Substituting (9) into Douglas' integral (2) and expanding in /3, we find that the linear 
term vanishes because 9*(a) is the minimum, while the quadratic part reads 

07T J Jo 1 - COS((7*(0i) - <7*(0 2 )J 

The function (5(9) obeys 

13(9) > -<j,(0) (13) 

for the derivative of the reparametrizing function to be positive. This is always satisfied if 
j3 is small and smooth enough. 

In order to calculate the semiclassical correction to the area law, we need to do the 
Gaussian integral 

[vp(9)e- Am0)] (14) 

with A 2 \p(0)\ given by Eq. (12). 

The typical values of f3, which are essential in the path integral over f3 in (14), are 
f3 ~ 1 / \J K S'min, i.e. small for V A'S'min ^> 1. Hence, the higher terms of an expansion 
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of A[9*(a) + (3] in (3 are suppressed [5] at large y/KS m - m - The loop expansion goes in the 
parameter l/KS min and only one loop contributes with the given accuracy. 

A comment is needed about the measure for the integration over (3(9). As is explained 
in Appendix A, the measure for the integration over a(8) involves a factor 

Because (3 ~ 1 / ^ K S niin , the second and third terms on the right-hand side are not essential 
to the given order. Therefore, the measure for the path integration over (5(9) in Eq. (14) is 
the usual one for smooth functions (3(9), while this factor will show up to the next order in 
1/ KS m i Q . 



III. PATH INTEGRAL OVER REPARAMETRIZATIONS: RECTANGLE 

In this section we show how the path integral over reparametrizations captures the 
Liischer term for a rectangle. 

The conformal map of the upper half plane onto the interior of a rectangle is given by 
the Schwarz-Christoffel mapping 

^b'")-' 2 — -■ < 16 > 

where 

F (*) = f «* (17) 

Vv^ / Jo V A* - x 2 a/1 - fix 2 
is the incomplete elliptic integral of the first kind. The two parameters A and \i are related 

to the coordinates of the vertices of the rectangle by 

AK(ii) = |, AK (yT=f) = T, (18) 

where K(jj) = F(l,/i) is the complete elliptic integral of the first kind. In deriving Eq. (18) 
we used the important identity 

F(l/n, fi) = K(jm) + iK (a/W^) • (19) 
Equations (18) relates fi to the ration of R/T as 



2T K 



R K(ji) 



; -. (20) 



In the limit T/R — > oo, when /i — > 0, this equation simplifies to 

T 4 

*5 = V (21) 

When z = s runs along the real axis, the variable u runs along the boundary of the 
rectangle with s = —l/^/JI,—y/jl,+y/JI,+l/y/jl (// < 1) mapped, respectively, onto the 
vertices of the rectangle: (-R/2,T/2), (-R/2, -T/2), (R/2, -T/2), (.R/2, T/2). The given 
choice of the argument of the mapping preserves the symmetry s — > 1/s. 

When z has positive imaginary part, the coordinates 

f z \ ( z \ AK 
X 1 (z) = AReF [^=^J , X 2 (z) = AlmF (^=,pj ^ (22) 

take their values inside the rectangle. These coordinates are conformal. For this reason we 
have 

/ s \ ( s \ AK 
x 1 {U{s))=AReF^—, f jLj, x 2 (U(s)) = AlmF ^— , fij ^ (23) 

whose implementation for the function t*(s) is discussed below. 

The boundary contour given by Eq. (23) satisfies Douglas' minimization (see Appendix B, 
Eq. (B3)). Correspondingly, the Douglas integral 

J —oo J — oo (,S & ) ^ ' 

as it should. We have verified these two equations numerically. 

A natural parametrization of the boundary of a rectangle is through r e S 1 : 

T T R R 

X\ — — tan r , x 2 = — arctan — < r < arctan — (25a) 

1 2 2 rp - - rp \ ) 

R R R R /nt-i \ 

Xi — — , x 2 = cot r arctan — < r < n — arctan — (25b) 

1 2 2 rp - - rp \ J 

T T R . . 

Xi — — tan r , x 2 = — n — arctan — < r < tt (25c) 

2 2 J. 

and analogously for negative r. Introducing 



t = tan-, (26) 
2 



we rewrite Eq. (25) as 



^i = J - , x 2 = < t < 27a 

1 — t 2 ' 2 R ~ ~ R v; 

Xl= «, l2 = R l^l ^f- f <K^ (27b) 



To relate t to s, we identify 



" :F(^=,A = T-^, (28a) 



tor - ^ < s < v> or < t < 

R f s dx T „t 2 -l 



. . - = R (28b) 

2K{n) ^^^T^ 2 At 

for yf\L < s < 1/y/fi or < t < . 

Solving the quadratic equation for t versus s, we obtain the minimizing function t*(s), which 
obviously obeys the boundary condition 

t(0) = 0. (29) 

The symmetry s — > 1/s plays apparently an important role. It guarantees that the points 
— oo, — 1, 0, +1, +oo are mapped onto themselves under the reparametrization s — > i*(s): 
**(-oo) = -oo, f*(-l) = -1, t*(0) = 0, = +1, t*(+oo) = +oo. 

It is convenient to invert Eq. (28) using the Jacobi elliptic functions. Inverting Eq. (28a), 
we get 

( 2K(fi)T t 
{ R Y^t 2 



^P^-^,/* , (30) 



_ ^ r- VT^TW - T JWTW - T 

for - V/i < s < V/i or < t < . 

-fx JrC 

The function sn has a nice trigonometric expansion in the parameter (nome) (Ref. [11], 
8.146.1) 



6XP \ n K(u) ) * W4) ^ h (31) 
giving 

s ~^Kf(rb))- (32) 

This formula is applicable for -R/2T < t < +R/2T, when —yfji < s < 
We can proceed in the same way with Eq. (28b), whose inverse is 

s = v^sn (k^i.) + iKfa) (^^ + |) ■ ( 33 ) 

Using the addition formula ([11], 8.156.1) and the reduction of sn(x, 0) and cn(a;, 0) to sin a; 
and cosx, we obtain 



s ~ ^//Icosh 



7T 



(34) 
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However, this expansion is useless for large t — > 2T/R, due to the imaginary part of the 
argument of sn the expansion will involve hyperbolic functions with arguments that can be 
large. Instead we can use the expansion of sn in terms of inverse sines ([11], 8.147.1), where 
these sines can be large, so only the first term is relevant: 

1 1 



, //; sin ( |+i^ ~ ~jf ) ) ^ ( '° Sh 



7T / i 2T\ 

-A*-* 1 -*) 



(35) 



Equation (34) is applicable for t — > R/2T from above and Eq. (35) is applicable when 
t ->■ 2T/R. 

The quadratic action, describing Gaussian fluctuations around t*(s), is like (12) rewritten 
for the real-axis parametrization: 

m m\ = S^ jM^k? VM - m)f • (36) 

The Gaussian approximations is justified for large KRT, when 

m ~ TWw- (37) 

To calculate the path integral over the quantum fluctuations around the minimizing 
function £*(s), we need a mode expansion of the infinitesimal reparametrizing function /3(t). 
To get rid of the projective symmetry, we keep fixed 3 points, e.g. —1, 0, +1 or —y/Jl, 0, y/Ji 
fixed: /3(-l) = 0, 0(0) = 0, = or 0(-^(-0J)) = 0, /3(0) = 0, P(U(y/p)) = 0. 

For each segment from ti to tf, we consider the mode expansion 



£(t) = ^c„sin ^^-^j 



(38) 



obeying the boundary condition /3(tj) = /3(t/) = 0. This set of sines forms a complete basis 
on the given interval. Actually we shall need the mode expansion for 4 segments attached 
to t — t*(± v //I) ~ ±R/2T because the large contribution of l//x will appear in A 2 only for 
those. 

The appearance of the large factor l//x is seen already from Eqs. (30), (33) because 
in the denominator in Eq. (36) is proportional to yfji. We need, however, to show that this 
l//i is multiplied by a factor ~ 1. 

Let us analyze the contribution to (36) that comes from —R/2T < t-i,t 2 < +R/2T, i.e. 
from the bottom side of the rectangle. Introducing the variable 

2T 

y=— t, -Ky<l (39) 



and using Eq. (32), we write the contribution of this domain to A 2 as 

^ r dVl f Am 7T jmzmf x n„c (40, 



4tt^ J_! J_! (sin(7ryi/2) - sm(ny 2 /2)y (*R 
with a positive constant C. A similar contribution appears if t x — R/2T ~ t 2 — R/2T ~ 
R 2 /(2T) 2 as is prescribed by Eq. (34) for both t\ > R/2T and t 2 > R/2T. An analogous 
contribution (with possibly some powers of T/R oc ln(4//z)) emerges also when t\ < R/2T 
and t 2 > R/2T or t x > R/2T and t 2 < R/2T. 

Therefore, the path integral over the quantum fluctuations around the minimizing func- 
tion t*(s) gives, using the ^-function regularization and Eq. (21) (and disregarding the logs 
to the order in consideration) 

"""""(if)- < 41 > 

modes \V^/ 



where the 4th power is due to the four sets of modes 1 . This remarkably reproduces the 
Liischer term (6) in d = 26! 

There was a subtlety in the derivation - the appearance of a logarithmic divergence at 
the corners of the rectangle if $ (corners) ^ 0. It is seen from Eq. (40), where the region 
near y x — y 2 — 1 or y x — V2 — —1 (associated with ti = t 2 — > R/2T or ti — t 2 — > —R/2T ) 
produces the logarithmic divergence 

d yi r dy 2 ^ (±1) = /3 2 (±1) In i (42) 
J -1+8 J -i+8 (yi + y 2 ^2) o 

with the upper (lower) sign referring to y = 1 (y = — 1). The coefficient is nonvanishing 
if /3(±1) 0. Analogously, the integral is logarithmically divergent at the corner, when 
fi,f 2 > R/2T or ij < i?/2T, t 2 > R/2T and vise versa. 

The logarithmic divergence can be regularized by smoothing the corners like in 
Refs. [7, 14]. It is clear from such a regularization that the contribution of trajectories 
with P(R/2T) 7^ to the path integral over /3(t) will be suppressed as the smoothing is 
removed. Consequently, this corner divergence does not effect the result of this section. In 
the next section we repeat the consideration for the case of an ellipse, when there are no 
corners. 

If d ^ 26, the asymptotic ansatz for the Wilson loops has to be improved to get the 
correct factor (d — 2)/24 in the Liischer term (6). This issue will be described in Sect. V. 



1 This is like in the computation of the static potential for the Polyakov string in Ref. [12]. 
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IV. PATH INTEGRAL OVER REPARAMETRIZATIONS: ELLIPSE 



In this section we evaluate the path integral over reparametrizations for an ellipse and 
obtain a prediction for the associated Luscher term. 

The necessary formulas are given in Appendix B of [6] and are partially reproduced in 
Appendix B. We are interested in the case of a very long ellipse when the ratio b/a — > and 
v — > 1 according to 

h i±> !*£gE^. (43 ) 



Using the asymptote 
we simplify Eq. (43) to 



2K{y) 
8 

m , . 

2 (l-i/)' 



K{v) v =$ Jin— ^— , (44) 



7T 2 



a 4 In 



(45) 



(l-i/) 

For v — > 1 the elliptic function simplifies and we have 



In 



2q+ v /4q 2 + (l- i /) 2 



fl.fr) = tt | ~ _|_ + 1| (46) 

for — 7r/2 < o" < 7r/2. Inverting Eq. (46), we find 

a,(0) = ^8hih(-ln JLV (47) 

Z \ 7T 1 — V J 

This is quite similar to Eqs. (32) and (34) for a rectangle with yjji replaced by (1 — v). 

The calculation of the path integral over reparametrizations at one loop is quite analogous 
to that for the rectangle described in the previous section. We see that the large factor of 
(1 — v)~ 2 emerges in Eq. (12) because cr* oc — from Eq. (47). To evaluate the coefficient, 
let us consider the domain of small fli and 9 2 : 



1 ,8 2 « [ ln—) , (48) 
which contributes to the integral in (12) 

modulo the powers of b/a oc ^ln (j~y ) • The same contribution comes also from the 
domain of both 9\ and fl 2 near n. We thus have 

A -- K jrhf (50) 
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for every mode. 

Integrating the Gaussian integral for every mode and using the ^-function regularization, 
we get (disregarding the logs to this order) a pre-factor of the type 

modes N ' N ' 

where the product runs over 4 sets of modes, which results in a factor of 4 in the exponent. 
This coincides with the Liischer term (6) for a rectangle of the size RxT in d = 26 dimensions 
provided 2 

- - - (52) 
R~ 2b- {b2) 

In Appendices C and D we confirm this by an explicit calculation of the determinant of the 
Laplace operator and the conformal anomaly for an outstretched ellipse. 



V. A GENERALIZATION TO ARBITRARY DIMENSIONS 

The results of two previous sections demonstrate the already mentioned fact that the 
ansatz (1) has to be modified in order to describe the Liischer term in d = 4 dimensions. 
A simple modification is based on the form of the path integral for a rectangle 

Jv(3{t)z- M[m]T o? ^ T/R (53) 

with quadratic action A 2 [/?(£)] given by Eq. (36). 

For an arbitrary curve this path integral can be expressed through the determinant of 
the corresponding operator, that enters A 2 , which we denote as O: 

J V(5(t) e - A *WW = (det 0)~ 1/2 . (54) 

It is now clear that the following modification of the ansatz (1) will provide the correct value 
of the Liischer term (6) in d dimensions: 

J Vt e ~ KA ^ (det 0)- (d - w)/4S . (55) 

This modification of Eq. (1) does not effect the classical limit, leading to the area law, and 
has the meaning of altering a pre-exponential in the semiclassical approximation. 

2 It is worth noting that the Wilson loops for a rectangle and ellipse then coincide [13] in d = 4 (the only 
dimension with &T/R Coulomb term) to the second order of perturbation theory. 
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To make the structure of the operator O more explicit, it is convenient to use the ex- 
pansion (9) of the direct function 9(a) rather than that (10) of the inverse function a(9) as 
above. Using the identity 

fi(U(s)) = --^fi( 8 ) 1 (56) 

ds 

which stems from the definitions (9) and (10), we then obtain for the real-axis parametriza- 
tion: 

m w)\ =£/<•«>/ d, 2 ^<£«M VM - j W r , (57) 

which determines the "momentum" (with respect to s) space operator 

0(pi,P2) = ^ J dq\q\(2x(p 1 +q)-x(-p 2 -q)-x(p 1 -p 2 +q)-x(-q)-x(q)-x(p 1 -p 2 -q) ] j (58) 
with 

x(p) = J dse ips x(U(s)). (59) 

We can finally substitute t*(s) by t(s) in this formula without changing the semiclassical 
approximation. 

It is worth noting that in contrast to the Laplace operator of Ref. [14], where the Liischer 
term was obtained for the Polyakov string, the present operator O lives in the boundary, 
which makes the construction nontrivial. 



VI. CONCLUSIONS 

The conclusion is that the reparametrization of the boundary curve involved in Eq. (1) 
carries information on the transverse fluctuations in 26 dimensions. As is shown in the 
previous section, it is possible to generalize this to any dimensions. 

Our motivation for the present paper is our previous work on the QCD/string scattering 
amplitudes [15], where we used that the amplitude can be expressed in terms of a Wilson 
loop through Feynman path integration. There we only considered the leading area behavior. 
However, having developed a path integral expression for the next term, we hope that the x^ 1 
integrals can be performed, thereby providing a momentum space analogue of the Liischer 
term. We hope this may help to answer a very interesting question as to how the intercept of 
the Regge trajectory changes under such a modification of the ansatz for the Wilson loops. 
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Appendix A: Measure for integrating over reparametrization 

Introducing 

Vi = Si-Si-i s N = s f , (Al) 
we rewrite the measure of [6] for the integration over reparametrizations as 

*"f ,. j-r f Si+1 dSi f S2 dSi 



/ ^difrs = lim | T f - — — — r / 

J So N ^°° to J SO ( S i+1 - S i) J SO 



'so j_2 J so V "H-i "%) J so 



S2 ~ - S ) 



i=l " /u 1 j=l 

The integration over v^s in Eq. (A2) can be represented through the integration over a 
scalar field as follows. Writing 

v i =e* i ' 2 , (A3) 

we have 

r°° c \ v . i r+°° 
Jo Vi 2 J_ 00 

and 

/ D diffS = Km _ J] / d^5^(s f -s -J2^ j/2 )- (A5) 

./so 00 z i=1 ./-oc J=1 

This represents the continuous measure as 

f ' V dm s = [ Vi, 6W (s f - So - dt e ~^ )/2 ) , (A6) 

J S() V ^ J SO 

where t is a certain parametrization of the contour (e.g. through the proper time) and Dip 
is the usual measure 

/ S f r+co 
= n / d ^( s ) • ( a? ) 

s=s J -°° 

The scalar field ip, that appears in Eqs. (Al), (A3), is in fact a discretization of the 
boundary value of the Liouville field 

<p(T,s) =i>(s), (A8) 

boundary 
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which is related to the boundary metric: 

^&y = „«.,/». (A9) 

Under the reparametrization s — > f(s) (df/ds > 0), the boundary value of the Liouville 
field ip(s) transforms as 

Ms) _> Mf( s )) = ^s) - 2 In ^ , (A10) 

that clarifies its relation to reparametrizations. 

The results of this Appendix make it possible to relate the ansatz (1) with Eq. (17) of 
Ref. [14], where the path integral over ifj is the same as the path integral over reparametriza- 
tion in view of Eq. (11) of [14]. They coincide in d = 26 provided S c \ is the Douglas integral. 



Appendix B: Douglas' algorithm for plane contours: conformal map 

The construction of the coordinates of the minimal surface, enclosed by a plane contour, 
is given by conformal mappings. Let us describe such a contour by two functions x\(t) and 
x 2 (t). Motivated by Appendix H of Ref. [4], we define the analytic functions 

Mlj J-oc vr {s'-zY ■ [t5L) 

The real and imaginary parts of $ M are 

f +OD ds' 

Re^(s) = / x^Uis')), (B2a) 



oo S ' ~ S 



» / v dxJtJs)) . , , .. dt*(s) . 

Im^(s) = M d * 1 = ^(^(s))-^. (B2b) 

Therefore, Douglas' minimization equation 

f + °° dgl ±(U(s)) ■ [x(U(s)) - x(U(sf))] = (B3) 

is satisfied if 

J]lm$J( S ) = (B4) 



at the real axis. 

For the circle we have 



Xl (t) + ix 2 (t)=C(t) = — (B5) 

1 ~~ r t 
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so Eq. (Bl) gives 



with 



$i(z) = i$ 2 (z) 



0, 



which obviously obeys Eq. (B4). ^From Eqs. (B6) and (B5) we find 

2s 



Im<E>i(s) = 
Im$ 2 (s) 



(1 + s 2 ) 2 ' 
1-s 2 



1 _ 2t 

2 Xl[t) ~ ~(TTW 
i . , , i-t 2 
2 X2[t) ~ (TT^F 



(1 + s 2 ) 2 ' 

and from Eq. (B2b) conclude that t*(s) = s, as it should be for the circle. 
For the functions obeying Eq. (B7), we always have <&i(s) = i$2(s) and 

$( s ) = $ 1 ( s )+i$ 2 ( s ) = 2$i(s); 

Im$i(s) + ilm$ 2 (s) = -^<3>(s). 

Equation (B2b) can then be rewritten as 

2 as as 



(B6) 
(B7) 

(B8a) 
(B8b) 



(B9a) 
(B9b) 

(BIO) 



For the given analytic function C(z) it determines the reparametrizing function t*(s). 
For the circle and ellipse it is more convenient to use the unit-disk parametrization, when 

du[c(u)-c(u>)\ (B11) 



= (u + If j 



27ri (u - u') 2 

For the ellipse C(oo) is given by the conformal map (where v is the same as s in Eq. (B26) 
of [6]) 



C(u) = V a 2 — b 2 sin 
and from Eq. (Bll) we obtain 



7T ^ ( u 
F -=,u 



$H = Va 2 - 6 2 



7T 



2AT(i/) Vv^ 7 ' 
(w + 1) 2 



(B12) 



2AT(z/) - wVl - z^ 2 



cos 



77 W" v 



(c 2 + 1) 



dcy) 

da; 



The final step is to extract d0*(cr)/dcr from Eq. (BIO). Remembering that 

C(6) = a cos 6 + i&sinfl, 

we find 

d6 1 



du - uj 2 \/\ - vu 2 



(B13) 

(B14) 
(B15) 



which for cu = e 1<T reproduces Eq. (B32) of [6]. 
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Appendix C: Calculating the Liischer term for an ellipse 



1. Elliptical coordinates and Mathieu functions 

Let us consider an ellipse as the boundary contour. To calculate the determinant of the 
Laplace operator with this boundary conditions, we parametrize the surface by elliptical 
(Lame) coordinates 

x\ = h cosh r cos a, X2 = h sinh r sin a" (CI) 

or 

X\ + 1x2 = h cosh (r + ia) . (C2) 

Here a e [0, 2ir) and r G [0, r ]. The boundary is approached for r = r = arctanh b/a and 
foci are at (±h, 0) with 



h = Va 2 - b 2 . (C3) 
The elliptic coordinates are conformal: 

ds 2 = h 2 (sinh 2 r + sin 2 a) (dr 2 + da 2 ) , (C4) 



so the area in the conformal gauge reads 

"to /•2ir 



(C5) 



The Laplacian reads in elliptic coordinates as 



1 f d 2 d 2 \ 

A= h 2 (sinh 2 r + sin 2 < 7) V^ 2+ ^ 2 J' ^ 

The Laplace equation separates in elliptic coordinates to the Mathieu equations 

d 2 F 

— 7 + (a-2gcos2cr)F = 0, (C7a) 
da 2 

d 2 C 

— -(a-2q cosh 2r) G = . (C7b) 
dr 2 

where a is a separation constant and q is related to the eigenvalue A by 

q = . (C8) 

A complete set of solutions for F and G is given, respectively, by the Mathieu functions 
ce m (a,q), se m (a,q) and the modified Mathieu functions Ce m (r, q), Se m (r,q) of integral 
order m. 
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FIG. 1: Plot of 1 — a m (q)/2q versus c = q/m 2 . 



While the characteristic numbers a m (q) and (3 m (q) for the ce m and se m modes are not 
explicitly known for large q ~ m 2 , which would be the case for an outstretched ellipse with 
6/a < 1, two asymptotic formulas exist [16] 



&m Pn 



mr^q 2 



1 q 2 5 q 4 9 q 6 1469 g 8 4471 q 



10 



m H — + 



+ 



+ 



+ 



2 m 2 32 m 6 64 m 10 8192 m 14 16384 m 18 



+ . 



(C9a) 



m 2 -€.q 1/2 1 2 1 m3 ^ m 4 33 m 5 63 m 6 527 m 7 

(C9b) 

These expansions are not formally applicable for c = q/m 2 ~ l, 3 but the first formula works 
numerically for c < 0.7, while second formula works numerically for c > 0.4 — 0.5, so there 
is an overlap. Special values of c are c = 0.61 and c = 0.64, when 



a m (cm 2 ) = (3 m (cm 2 ) = 2cm 2 = 2q . 



(CIO) 



They are obtained numerically with Mathematica, but are derivable from the expansions 
(C9a) and (C9b), correspondingly. As we see from Fig. 1, the difference 2q — a m (q) is positive 
for c > 0.64, where the expansion (C9b) is applicable. 

The boundary condition requires Ce m (r ) = and Se m (r Q ) = 0, which implies 



T — %r. 



(Cll) 



where Z m ^ n are the corresponding zeros. This determines the eigenvalues of the Laplacian 
with m and n playing the role of angular and radial quantum numbers, respectively. For 



Some results concerning the Mathieu characteristic numbers for large q ~ m 2 can be found in Ref. [17]. 
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small r (= large a/6) we have [16] 



so that 



Ce m (r) Rj const, x cos (^y / 2q — a m T^ (C12a) 
S'e m (r) « const, x sin (a/ 2 <? - /3 m rj (C12b) 



v/2g - a m r = — (C13) 



because a m = (3 m for large m. Here n is odd or even for the Ce or Se modes, respectively. 
We conclude therefore that 

/7ina\ 2 

2q ~ am= V2b) (C14) 
has to be large and justify large q ~ m 2 . 

2. Evaluating the determinant 

Each mode contributes 

det (A- 1 / 2 ) = l[\-X 2 = e E "- ln ( A "-) . (C15) 

n,m 

Since a is large, the sum over m can be replaced by an integral over oj = m/a like in Ref. [18] 
and for large q 1 / 2 ~ m ~ a/6, we have 

K!l = rf Q , (C16) 

where 



so that 



£ ln ( A i'. 2 ) = -<■ /°° d ^ k [ r f (7)] • < C18 > 



For large r the integral on the right-hand side is proportional to r and we can get the 
coefficient of proportionality by differentiating with respect to r. This gives 



POD 

^ (K!l 2 ) = -ar / dx [1 - xf(x)/f(x)] 
Jo 



(C19) 
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1 .95 



2.05 



FIG. 2: Matching the expansions (C23) and (C24). The region near x = 1.95 is enlarged in the 
right figure. 



and finally we obtain 

£ ln ( A ^ 2 ) = /"^[l -*/'(*)//(*)] 

n,m n ^ 

= Y^h\ dx i 1 - x f'( x )/f( x )]- 

J 



(C20) 



For an T x R rectangle with T R, when 



A, 



7i 2 n 2 7r 2 m 2 



+ 



R 2 T 2 



rectangle 



resulting in 



f(x) = yl + x 2 rectangle 



(C21) 



(C22) 



we reproduce the well-known result for the Luscher term by the substitution a = T/tt, 
b = R/2. 

For an outstretched ellipse we have from the expansion (C9b) 



f I = 1+xH h 

J v 7 4 16 128 1024 2048 

for x < Xo ~ 1.95, where the series converges, and 

r/ x *>*o , ^ „ _ 0.960520 1.041899 2.418227 7.138640 23.73152 
j(x) = 1.600095a; H 1 ; - 1 - 



(C23) 



x 



X' 



X ' 



X a 



84.70585 



x 



ii 



+ . . . 



(C24) 



for x > Xq ~ 1.95. They match each other pretty well, as is depicted in Fig. 2. The numerical 
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value of the integral on the right-hand side of Eq. (C20) is then 2.84 to be compared with 
7r as is expected from Eq. (51). 

To clarify the calculation, let us keep only three terms on the right-hand side of Eq. (C9b), 
when we get a quadratic equation and analytical formulas are available. We then find 



f(x) =x+ y/l+ x 2 /2 (C25) 

and expanding it either in x or in l/x we obtain analogs of Eqs. (C23) and (C24). The 
two expansions now match at xq = V2 = 1.41 to be compared with x = 1.95 for all 
nine terms left on the right-hand side of Eq. (C9b). The value of the integral is now 
a/2 log [3 + 2^2] =2.49 which is smaller than 2.84, which may characterize to which extend 
the approximation of a m by the eight term in (C9b) is better than by the tree terms. In 
Appendix D we confirm the value of tt by another method. 

Appendix D: Liischer term from conformal anomaly 

A starting point is the formula [7, 19] 

Liischer term = / drdcr d a In p d a In p (Dl) 

96tt J a y a y v J 



that relates the Liischer term to the conformal anomaly. Here the metric p is 

u = t + icr (D2) 



P = 
where 



dz 2 



du 



co = h cosh (r + icr) (D3) 

runs inside the ellipse with r and a being elliptic coordinates as in Eq. (C2). The function 
(inverse to (B23) from [6]) 

, , r- (IK (tx \\ r- / 2K , A r- cn [u - it)) 

z(u) = v^sn ^— (- - a + ir ) ) = v^sn (tf - — (a - it) ) = ^ - {2J£ ^ _ ^ 

(D4) 

conformally maps ellipse onto a unit disk. 

To calculate (D2) for a/b 3> 1, we substitute 

cn(Mi u>I/ ) _ i 



dn(^m,z/) dn(^w,z/') 
20 



(D5) 



with 

v' = Vl - v 2 . (D6) 

Differentiating, we get 

dw vr dn 2 (^u, i/) 

Using the fact that Jacobi elliptic functions reduce to trigonometric functions as v' — > and 
substituting 

K^p y (D8) 



we infer from Eq. (D7) 



and 



dz ( . Ti a \ ( a , 

— oc exp — i— -u I ==>- p oc exp 7r-cr I (D9) 
du V 2 o / ' V 6 ■ 



Liischer term = / dr da [tt-] = — - (D10) 

967T70 J \ bJ 48 6 V ; 



that confirms the extra it/ 2 in Eq. (51) for an ellipse in comparison with a rectangle. 
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